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1.  A PUZZLING  PATTERN 

Following  the  development  and  installation  of 
a Pulse  Pair  Processor  (Novick  and  Glover,  1975) 
and  a color  display  (Jagodnik  et  al.,  1975)  for 
analysis  and  presentation  of  Doppler  spectral 
moments,  we  have  routinely  acquired  velocity  pat- 
terns in  all  available  storms,  using  the  5.5-ca 
Porcupine  Doppler  radar.  The  usual  procedure  for 
taking  data  requires  the  antenna  to  scan  in  azimuth 
at  a constant  elevation  angle.  After  each  complete 
rotation,  the  antenna  elevartion  is  stepped  up  to 
a new  value.  The  elevation  increments  are  general- 
ly 1°,  commensurate  with  the  antenna  half-power 
beamwidth  of  0.9°,  hut  often  larger  Increments  are 
used  above  an  elevation  of  10°.  The  output  of  the 
Pulse  Pair  Processor,  containing  information  on 
reflectivity,  velocity  mean,  and  velocity  variance, 
in  each  resolution  cell,  is  recorded  on  tape  for 
subsequent  analysis  and  is  presented  in  real  time 
on  the  color  display  in  FPI- format.  The  color 
display  of  data  taken  at  any  elevation  angle 
above  0°  is  not,  strictly  speaking,  a plan  view  of 
the  storm.  It  is,  instead,  a conical  surface  with 
apex  at  the  radar,  slicing  across  the  storm  at 
heights  defined  by  range  and  elevation  angle. 

The  summer  of  1975  in  eastern  Massachusetts 
was  notable  for  its  scarcity  of  severe  thunder- 
storms and  for  occasional  intrusions  of  widespread, 
winter-type  precipitation.  A small  coastal  north- 
easter, a weak  cousin  of  the  producers  of  heavy 
vinter-time  snow  along  the  Atlantic  seaboard, 
visited  our  area  briefly  during  the  morning  of 
August  7.  Precipitation  was  light,  with  an  average 
rate  of  only  0.6  mm/hr,  but  was  sufficient  to 
provide  detectable  tracers  of  velocity  for  our 
Doppler  radar.  The  storm  vas  scanned  in  ^-incre- 
ments of  elevation  angle  from  1®  to  10°. 

A photograph  of  the  velocity  pattern  at  eleva- 
tion 8°,  displayed  in  color,  is  reproduced  as 
Figure  D of  the  color  plate  in  the  preceding  paper 
by  Kraus  and  Donaldson.  As  in  the  other  figures  of 
this  plate,  the  colors  represent  contours  of 
velocity,  in  accordance  with  the  scale  on  the  right 
of  the  photograph,  ranging  from  ♦22  m/s  away  from 
the  radar  (red),  down  through  yellows  and  greens 
to  black  (±  1.3  m/s)  and  on  to  negative  velocities, 
toward  the  radar,  represented  in  succession  by 
white,  ever-deeper  shades  of  blue,  and  finally 
purple  indicating  velocities  up  to  -22  m/s.  The 
narrow  black  band  separating  advancing  and  receding 
velocities  not  only  eliminates  the  stationary 
ground  clutter,  but  also  serves  as  an  indicator  of 
locations  where  motions  within  the  storm  are  normal 
to  the  radar  beam.  The  two  range  rings  are  located 
at  16  and  32  km. 

The  pussling  feature  of  Figure  D,  and  its  neigh- 
bors at  nearby  elevation  angles,  is  the  asynoetry 


of  the  color  patterns.  Maximum  positive  and  negative 
velocities  are  sensed  by  the  Doppler  radar  when  its 
antenna  beam  is  pointing  downwind  and  upwind.  In  a 
horizontally  homogeneous  wind  field,  archetypical  of 
widespread  stratiform  storms,  we  would  expect 
Doppler-sensed  velocities  to  decrease  uniformly  on 
either  side  of  a maximum,  in  accordance  with  the 
cosine  of  the  directional  change  of  the  beam  from 
the  velocity  maximum.  Curiously,  though.  Figure  D 
shows  a concentration  of  velocity  contours  near 
each  peak,  but  each  has  a long  tail  trailing  clock- 
wise. The  pattern  resembles  the  Yin-Yang  symbol  of 
Chinese  cosmology,  or,  in  a less  traditional  orienta- 
tion, two  tadpoles  chasing  one  another  in  a circle. 


2.  SOLUTION  OF  THE  PUZZLE  BY  FOURIER  ANALYSIS 

No  ravinsonde  winds  were  available  on  August  7, 
1975  for  comparison  vith  the  Yin-Yang  Doppler  velo- 
city pattern.  Consequently,  we  generated  a wind 
field  from  the  recorded  velocity  data,  using  the 
precipitation  as  tracers  and  following  the  scheme 
for  Fourier  analysis  of  harmonics  developed  by 
Browning  and  Vexler  (1968).  Our  aim  was  not  only 
to  derive  a vertical  profile  of  average  wind,  but 
also  to  gain  some  knowledge  of  the  uniformity  of  the 
wind  at  each  height.  At  the  time  of  this  writing, 
the  harmonic  analysis  is  confined  to  a quasi- 
cylindrical  section  through  the  storm  with  a radius 
of  8.325  km  around  the  radar.  This  section  extends 
up  to  a height  of  1.1*5  k».  using  integral  values  of 
elevation  angle  from  1°  to  10°. 


Browning  and  Vexler  suggested  that  the  Fourier 
coefficients  of  order  0,  1,  and  2,  which  give  the 
horizontal  wind  field  properties  of  divergence,  mean 
wind  velocity,  and  deformation,  respectively,  may  be 
computed  with  sufficient  accuracy  by  summations  of 
Doppler  velocities  at  10°  intervals  of  azimuth. 
Accordingly,  we  sampled  the  velocities  at  36  equally- 
spaced  azimuths  around  the  scanning  circle,  for  each 
of  the  eight  elevation  angles  from  2°  through  9°- 
At  elevations  1°  and  10°  ve  could  obtain  data  at 
only  30°  azimuth  intervals,  owing  to  a data  gap  at 
elevation  10°  and  several  places  at  1°  where  ground 
clutter  strongly  contaminated  the  precipitation  echo. 


At  these  two  extremes  of  the  elevation  angles  we 
computed  only  the  mean  wind  vector,  because  we  felt 
that  a 30°  azimuth  spacing  was  much  too  coarse  for 
meaningful  estimates  of  divergence  and  deformation. 

Ve  will  not  lay  out  the  computational  details  here; 
the  reader  who  may  be  interested  in  performing 
harmonic  analyses  of  Doppler  velocities  will  find  all 
the  instructions  he  needs  in  the  treatise  of  Browning 
and  Vexler. 

Estimates  of  divergence  using  Doppler  velocities 
acquired  at  elevation  angles  above  0°  are  contamina- 
ted by  a small  component  of  precipitation  fall  speed 
moving  toward  the  radar.  Ve  did  not  measure  the 
vertical  speed  of  the  raindrops,  but  on  the  basis  of 
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the  rain  rate  (0.6  am/hr)  ve  estimated  a minimum 
value  of  3 m/s  for  fall  speed.  Since  all  diver- 
gence values  were  negative  by  an  amount  greater 
than  the  fall  speed  contamination,  ve  feel  ve  have 
made  maximum  estimates  of  convergence.  Although 
this  procedure  lacks  the  rigor  requisite  for  a 
dynamic  analysis  of  the  vind  field,  ve  feel  it  is 
more  than  sufficient  for  our  purpose  in  assessing 
the  uniformity  of  the  vind. 


Table  I 

Wind  Field  Properties  of  7 August  1975  Storm 
Computed  by  Fourier  Analysis 


■leva-  Wind  Direo-  Diver-  Defor- 
tion  Hgt.  Speed  tion  gence  nation 
mgle  (m)  (m/s)  From  (10-5  a-1) 


Results  of  the  vind  field  analysis  are  listed 
in  Table  I,  and  vind  vectors  are  plotted  in  hodo- 
graph  form  in  Fig.  1.  We  nov  have  the  ansver  to 
tvo  questions:  What  is  the  vertical  profile  of 

mean  vind,  and  hov  homogeneous  is  it  at  any  given 
height? 


up  to  the  speed  maximum  at  870  m,  the  vind  veered 
about  12°,  averaging  a direction  change  of 
2.1°/100  m.  Above  the  speed  maxi mum  the  vind 
veered  much  more  sharply,  at  a rate  of  7.3°/100  m, 
until  at  1450  m it  vas  bloving  directly  from  the 
east.  At  higher  altitudes  the  vind  eventually 
became  vesterly,  indicated  by  the  green  (outbound) 
patch  in  Fig.  D tovard  the  east  betveen  the  16  and 
32  km  range  rings.  Examination  of  the  recorded 
velocities  directly  east  of  the  radar  at  elevation 
snglss  7°  through  10°  shoved  velocities  changing 
from  Inbound  to  outbound  at  a height  of  2.7  to 
2.8  km.  This  altitude  caps  the  easterly  regime  and 
is  the  base  of  the  overlying  vesterlles. 

We  nov  have  all  the  information  ve  need  to 
decipher  the  kinematics,  if  not  the  cossiology,  re- 
vealed in  the  color  display  of  velocities.  At  any 
elevation  angle  a vfaere  the  radar  can  locate  targets 
above  the  vindapeed  maximum,  there  vil 1 be  a locus, 
or  spine,  of  maximum  velocities  in  both  the  out- 
bound and  inbound  halves  of  the  picture.  The 
absolute  maximum,  or  if  you  vill,  the  peak  of  the 
spine,  occurs  at  a particular  range,  rB,  vhere 
rB  sin  a ■ 870  m,  and  at  azimuth  angles  tovard  the 
northeast  (inbound)  or  northwest  (outbound),  vhich 
are  directions  vhere  the  radar  beam  lies  along  the 
peak  vind.  Clockwise  from  the  peaks,  the  radar 
encounters  winds  vhich  are  moving  along  the  same 
azimuth  as  the  beam  at  a greater  height  (as  indi- 
cated in  Fig.  1)  and  hence  at  a greater  range. 

Since  vind  speeds  are  decreasing  slightly  in  this 
altitude  interval,  the  indicated  waxinnm  will  be  at 
a compromise  range  (height)  somewhere  betveen  the 
vind  direction  at  that  height  and  the  somewhat 
greater  speeds  found  at  a lower  height.  The  spine 
of  maximum  velocities  consequently  spirals  outward, 
tovard  higher  ranges,  vith  Increasing  azimuthal 
distance  clockwise  from  the  absolute  peak  velocity. 

On  the  other  side,  i.e.,  counter-clockwise  from 
the  absolute  peak  velocity,  the  same  general  struc- 
ture is  found  but  in  a strikingly  asjnuetrical  manner 
In  this  region  the  radar  beam  is  directed  along  winds 
at  heights  below  the  absolute  velocity  peak,  so  the 
velocity  spine  tends  to  continue  its  Inward  spiral 
vith  counterclockwise  direction.  However,  the  vind 
field  is  radically  different  here:  vith  a <mii 
decrease  of  azimuth  the  height  of  maximum  vind  in 
that  direction  and  its  maximum  speed  both  decrease 
precipitously.  Accordingly,  the  velocity  contours 
in  this  region  are  crowded. 


Fig.  1.  Hodograph  representation  of  vind  field  of 
7 August  1975  computed  by  Fourier  analysis  of 
Doppler  velocity  components.  Numbers  adjacent  to 
every  other  dot  shov  heights  in  meters  above  ground 
level.  The  surface  datum  vas  recorded  by  anemo- 
meter. 


Considering  the  second  question  first,  ve  see 
that  values  of  divergence  and  deformation  are  in 
the  vicinity  of  10" 5 to  1(TU  s”1,  typical  of  wide- 
spread precipitation  occurring  in  synoptic-scale 
disturbances,  and  at  least  tvo  orders  of  magnitude 
less  than  shear  values  associated  with  severe 
thunderstorms.  The  average  value  of  deformation, 
multiplied  by  the  radius  of  the  area  for  which  it 
was  computed,  is  only  7-4J  of  the  average  wind 
speed.  For  divergence,  the  comparable  value  is 
only  3.6jf  of  average  vind  speed.  We  are,  therefore, 
satisfied  that  ve  have  a vind  field  vhich  is  rather 
close  to  a condition  of  horizontal  homogeneity. 

For  the  purpose  of  explaining  the  arresting  features 
of  the  color  display,  ve  need  consider  only  the 
vertical  profile  of  mean  wind. 


The  vind  speed  Increased  from  5 m/s  at  the 
surface  (measured  by  anemoawter)  up  to  a maximum 
of  18.8  m/s  at  the  870  m level.  The  direction  of 
the  vind  vas  nearly  constant  from  slightly  north  of 
northeast  over  the  lowest  290  m.  From  that  level 


In  sum,  our  tvo  tadpoles  have  tails  because  the 


wind  profile  has  a marked  asymmetry  on  either  side 
of  its  speed  maximum.  On  the  "tail"  side  the 
speed  shear  is  small  hut  the  directional  shear  is 
large.  Note,  however,  that  the  tadpoles  are  not 
aligned  in  a circle,  which  would  he  out  of  charac- 
ter with  a horizontally  homogeneous  wind  field. 
Rather,  they  are  spiralling  inward.  Note  also 
that  at  any  given  range  the  measured  velocity  com- 
ponents are  symmetric  with  respect  to  azimuth  angle 
on  either  side  of  a velocity  maximum.  This  condi- 
tion is  required  hy  a uniform  wind  field  at  any 
height . 

3.  COLLATERAL  ANALYSES 

Contemplation  of  the  asymmetries  in  the 
Yin-Yang  velocity  pattern  played  a part  in  inspira- 
tion of  the  Kraus  and  Donaldson  (1976)  work  on 
generation  of  velocity  patterns  from  model  wind 
fields.  Figure  6 in  that  study,  shoving  the 
velocity  pattern  associated  with  veering  wind 
direction  and  speed  increasing  up  to  a maximum  and 
then  decreasing  above  it,  comes  closest  to  a depic- 
tion of  our  Fig.  D.  However,  our  pattern  is  much 
less  symsetric  because  of  the  strikingly  different 
character  of  our  vind  field  above  and  belov  the 
speed  maximum. 

Novick  (1976)  was  also  Intrigued  by  the  novel 
appearance  of  our  velocity  pattern.  He  focussed 
his  concern  on  the  nature  of  the  S-shaped  zero 
velocity  line,  vhich  appears  as  a black  band 
separating  the  inbound  and  outbound  colors  in 
Fig.  D.  He  derived  an  analytical  solution  for  the 
equation  of  this  line,  under  the  assumption  of 
constant  wind  speed  and  a linear  change  in  vind 
direction  vlth  height.  He  found  that  the  zero- 
velocity  curve  can  be  expressed  simply  as  a spiral 
vlth  distance  from  the  origin  linearly  proportional 
to  azimuth  angle,  nils  is  certainly  very  similar 
to  the  picture  we  see  at  heights  above  the  speed 
maximum,  where  our  computed  wind  speed  decreased 
by  only  22%  while  the  wind  direction  veered  1»2°. 
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In  the  course  of  our  study  another  thought 
occurred  to  us.  A complex  pattern  of  velocity  may 
be  Interesting  not  only  for  what  it  may  reveal, 
but  also  for  what  it  may  conceal.  In  a situation 
amenable  to  Fourier  analysis  of  the  Doppler  veloci- 
ties, or  one  where  environmental  winds  are  avail- 
able from  another  source,  we  could,  in  principle, 
generate  a display  which  would  present  the  observed 
Doppler  velocities  minus  the  horizontally  uniform 
vind  components  along  the  radar  beam.  Whatever 
remains  on  display  after  such  treatment  would  be 
solely  a function  of  anomalous  winds  indicative  of 
distortions  such  as  vortices.  Jets,  waves,  reac- 
tions to  obstacles,  and  the  like.  These  anomalies,' 
having  special  value  for  severe  storm  forecasting,^* 
would  stand  out  in  isolation  once  the  masking  V 
uniformities  are  removed.  ’ 
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